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■ Abstract 

We study the lepton number violating , e"*" ) conversion in nuclei mediated by the 
\^ ' exchange of virtual light Majorana neutrinos. We found that a previously overlooked 

imaginary part of this amplitude plays an important role. The numerical calculation 
has been made for the experimentally interesting ,e'^) conversion in ^^Ti using 
^ \ realistic renormalized proton-neutron QRPA wave functions. We also discuss the very 

' similar case of the neutrinoless double beta decay of ^®Ca. The ratio of (//"jC"*") 

(<~^ . conversion over the total /x~ absorption has been computed taking into account the 

\ current constraints from neutrino oscillation phenomenology. We compare our results 

■ with the experimental limits as well as with previous theoretical predictions. We have 
found that the Majorana neutrino mode of e"*") conversion in '^^Ti is too small to 

^ , be measurable in the foreseeable future. 

1 Introduction 

>■ . Lepton number (L) conservation is one of the most obscure sides of the standard model(SM) 
^ ' not supported by an underlying principle and following from an accidental interplay between 
^ \ gauge symmetry and field content. Any deviation from the SM structure may introduce 
■ - - ' L non-conservation(^). Over the years the possibility of lepton number non-conservation 
has been attracting a great deal of theoretical and experimental efforts since any positive 
experimental ^ signal would request physics beyond the SM. In addition it would also show, 
that neutrinos are Majorana particles 

Recent neutrino oscillation experiments practically established the presence of non-zero 
neutrino masses, a fact that itself points to physics beyond the SM. However neutrino oscil- 
lations are not sensitive to the nature of neutrino masses: they can either be Majorana or 
Dirac masses leading to the same observables. 

The principal question if neutrinos are Majorana or Dirac particles can be answered only 
by studying a lepton number violating processes since AL = 2 is a generic tag of Majorana 
neutrinos. Various lepton-number violating processes have been discussed in the literature in 
this respect (for recent review see 0). They offer the possibility of probing different entries of 
the Majorana neutrino mass matrix Mlj\ Among them there are a few ^ nuclear processes 
having prospects for experimental searches: neutrinoless double beta decay (Oi//5/5), muon 

^On leave of absence from the Joint Institute for Nuclear Research, Dubna, Russia 
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to positron (/i~,e+) conversion and, probably, muon to antimuon , ji'^) conversion 0. 
They probe M^^^^ M^^^ and matrix elements respectively. Currently the most sensitive 
experiment intended to distinguish the Majorona nature of neutrinos are those searching 
for neutrinoless Oi//3/3-decay [||, ||, ^ The nuclear theory side |] of this process have 
been significantly improved in the last decade allowing reliable extraction of fundamental 
particle physics parameters from experimental data. Muon to positron nuclear conversion 
(/i~,e"'") is another ^ nuclear process with good experimental prospects. 

The important role of muon as a test particle in the search for new physics beyond the 
standard model has been recognized long time ago. When negative muons penetrate into 
matter they can be trapped to atomic orbits. Then the bound muon can disappear either 
decaying into an electron and two neutrinos or being captured by the nucleus, i.e., due 
to ordinary muon capture. These two processes conserving both total lepton number and 
lepton flavors have been well studied both theoretically and experimentally. However, there 
are two other not yet observed channels of muon capture: The muon-electron (/i~,e~) and 
muon-positron e"^) conversions in nuclei 0, |T^, [T3|, |T5[: 



{A,Z)+iil ^ e- + {A,Z)\ (1) 
{A,Z)+ii-, e+ + (A,Z-2)*. 

Apparently, the (/i~,e^) and (/i~,e~) conversion processes violate lepton number L and 
lepton flavor Lf conservation respectively. Additional differences between the (/i~,e~) and 
(/x~,e^) lie on the nuclear physics side. The first process can proceed on one nucleon of 
the participating nucleus while the second process involves two nucleons as dictated by 
charge conservation [|1^, Note also that the (/x~,e~) conversion amplitude is quadratic 
and (/i~,e^) amplitude linear in the neutrino mass. Thus the second process looks more 
sensitive to the light neutrino masses. The present experimental limit on e"*") conversion 
branching ratio in ^®Ti is [|T6| , p!7| 

R,e^{Ti) = . ,..oJ < 4.3 X 10-^1 (2) 



In the present paper we study the light Majorana neutrino mechanism for the {fi~,e~^) 
conversion. Despite the previous rough estimates 0] indicate a very small branching ratio 
for this mode of e"*") conversion, by far below the experimental bound (^, we undertake 
a detailed study of this mode for several reasons. First, the nuclear theory of {fi~,e'^) 
conversion is not yet sufficiently elaborated, as in the case of Oi^/9/9-decay, and requires further 
development. Second, (/i~,e+) conversion may receive contribution from other mechanisms 
offered by various models beyond the SM such as the R-parity violating supersymmetric 
models, the leptoquark extensions of the SM etc. Some of these mechanisms may involve the 
light neutrino exchange and, therefore, from the view point of nuclear structure calculations 
they resemble the ordinary light neutrino mechanism. Thus our present study can be viewed 
as a first step towards a more general description of (/i~,e+) conversion including all the 
possible mechanisms. 

Below we develop a detailed nuclear structure theory for the light neutrino exchange 
mechanism of this process on the basis of the nuclear proton-neutron renormalized Quasi- 



particle Random Phase Approximation (pn-QRPA) wave functions [T^, |l^]. We perform a 
realistic calculation of the width of this process for the nuclear target ^^Ti using limits on 
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neutrino masses and mixings from neutrino oscillation phenomenology. A comparison with 
the previous estimations of -R^e+ will be also presented |Tl|. 

The paper is organized as follows. The possible values of Majorana neutrino masses and 
mixings are discussed in sect. 2. The amplitude and width of e"^) conversion are derived 
in sect. 3. The details of the calculation for the case of e"*") conversion in '^^Ti and our 
results are given in sect. 4. In sect. 5 we summarize our conclusions. 



2 Majorana neutrino mass matrix 

The finite masses of neutrinos are tightly related to the problem of lepton flavor/number 
violation. The Dirac, Majorona and Dirac- Majorana neutrino mass terms in the Lagrangian 
offer different neutrino mixing schemes and allow various lepton number/flavor violating 
processes pO| , pT| , p2| . The favored neutrino mixing schemes has to accommodate present 
neutrino phenomenology in a natural way, in particular, to answer the question of the small- 
ness of neutrino masses compared to the charged lepton ones. The most prominent guiding 
principle in this problem is the see-saw mechanism which can be realized in various models 
beyond the SM. A generic neutrino mass term is given by the formula 

^D+M ^ _ ^ [1 ^^M)^,^ + ^ ^ {M^% (^;«)^+^ {M%i u[,] + h.c. (3) 

The first two terms do not conserve the total lepton number L. Here, z/^ and z/^ are the 
weak doublet and singlet flavor eigenstates. The indices L and R refer to the left-handed 
and right-handed chirality states, respectively, and the superscript c refers to the operation 
of charge conjugation. and are complex non-diagonal symmetrical 3x3 matrices. 
The flavor neutrino fields are superpositions of six Majorana fields Vi with definite masses 
rrii. Yanagida, Gell-Mann, Ramond and Slansky suggested that the elements of and 
be comparable with the masses of charged leptons and the hypothetical scale of lepton 
number violation (Mj;,jvy ^ IQ^'^GeV), respectively. Then by diagonalization of the Dirac- 
Majorana mass term one ends up with the three very light and three very heavy neutrino 
eigenstates. This is the celebrated see-saw mechanism. Enlarging the number of the right 
handed neutrino states vr one can introduce sterile light mass eigenstates which may play a 
certain role in the explanation of the neutrino oscillation data including the LSND results. 
However the active-sterile neutrino oscillations as a dominant channel seems to be disfavored 
according to a recent Super Kamiokande global analysis [^] and work in Tiibingen |3^ . 



Sticking to the three neutrino scenario one may try to reconstruct the corresponding 
mass matrix from the neutrino oscillation data. This requires certain assumptions on its 
structure or additional experimental data. Solar, atmospheric and LSND neutrino data give 
information on the neutrino mass square differences Am?- as well as on the mixing angles 



of the unitary matrix U [0, ^ ^ ^ ^ relating the weak z/^^ and mass Vil neutrino 
eigenstates 



IL 



Y: l^^L {I = e, /X, r). (4) 



This information can be used to restore the neutrino mass matrix inverting its diagonalization 



as 



P[ph ^ . (iiag{mi, 1712, 1713,) ■ t/^'^)^, (5) 
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if additional assumptions about the overall mass scale as well as about the CP phases (^^p 
of the neutrino mass eigenstates are made. This matrix can be identified with the 3x3 
Majorana mass terms M^^ in (|^) and used in various phenomenological applications, for 
instance, in analysis of lepton number violating processes. The recent literature contains 



many sophisticated studies made in this direction (see, for instance, and references 
therein) . 

The elements of the Majorana mass matrix are related to the effective Majorana neutrino 
masses < >a/3 (a, /3 = e, /x, r) as 

M% ^< = j: Uil^ f/(? C^S m,, (6) 

k 

neglecting mixing with heavy neutral states if they exist in the neutrino mass spectrum. Am- 
plitudes of lepton number violating processes are proportional to the corresponding effective 
neutrino masses p|, Thus Oz//3/3-decay amplitude is proportional to < rrii, >ee, the so 



called effective electron neutrino mass [|^ From the currently most stringent lower limit 
on Oz//?/5-decay half-life of '^^Ge T^J^ > 1.1 x W^years one obtains < >ee< 0.62 eV 
The effective Majorana muon neutrino mass < irty is related with the light 



neutrino exchange modes of muonic analog of Oz//?/3-decay [Q, semipletonic decay of kaon 
— »• 7i~^^^~^ PH , etc. < rrii, >^e enters the amplitude of the {fi~,e~^) conversion 



lOl , [TT[| and of the kaon decay into a muon and a positron (K^ n fi^e~^). Some other 
elements of Ai^^ are associated with rare B-decays p3|. 




In Ref. the maximal allowed values for the elements of Ai^^ have been deduced from 
solar, atmospheric, LSND data and restriction coming from Oi//3/3-decay. The result is 

< >e/. < >er \ ( 0.60 0.97 0.85 \ 

< >MM <^u>^,r \ < 0.97 0.76 0.80 eV. (7) 

< < >rr J V 0.85 0.80 1.17 / 

Thus in our analysis of the light neutrino exchange mode of the e+) conversion process 
we shall assume | < rrii, >^e | < 0.97 eV. 

We shall also use a more conservative model independent estimate of the effective neutrino 
mass. Atmospheric and solar neutrino oscillation data show up Am^ << (leV)"^ suggesting 
that all the neutrino mass eigenstates are approximately degenerate at the 1 eV scale |PS 



This observation in combination with the tritium beta decay endpoint allows one to set upper 
bounds on masses of all the three neutrinos Tne^^^r < 3eV. Thus in the three neutrino 



scenario one derives | {mu)ij \ < 9eV for i,j = e, /x, r P, |3^ . 



3 The (/i , e+) conversion mediated by light neutrinos 

The process of (/i^, e+) conversion is very similar to the Oi^/?/?-decay. Both processes violate 
lepton number by two units and take place only if the neutrino is a Majorana particle with 
non-zero mass. However, there are other important differences: i) The available energies 
for these two processes differ considerably. In addition, the number of leptons in final 
states is different. These facts result in significantly different phase space integrals, ii) The 
emitted positron in (/i~,e+) conversion has large momentum and therefore the long-wave 
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Figure 1: The direct (a) and cross (b) Feynmann diagrams of the (/i ,e~^) conversion in 
nuclei mediated by Majorana neutrinos. 



approximation is not vahd. iii) As it will be shown below, in the case of light neutrino- 
exchange there is a singular behavior of the (/i~,e~*") nuclear matrix element which is the 
additional source of the difficulties for the numerical integration, iv) In the case of the 
(/i~,e+) conversion there is a great number of nuclear final states. Nevertheless, the major 
contribution is here assumed to come from the transition to the ground state of the final 
nucleus. 

We shall discuss the amplitude and width of (/i~,e"'") conversion in nuclei mediated 
by light Majorana neutrinos. This process is shown in Fig. |1|. We concentrate only on 
the nuclear transition connecting the ground states of the initial and final nuclei, which is 
favored from the experimental point of view due to the minimal background. In this case 
the e"*" spectrum has a peak at the energy 

E,+ =m^-eb-{Ef-E,). (8) 

Here, m^, Sb, Ei and Ef are the mass of muon, the muon atomic binding energy (for '^^Ti 
Eb = 1.45 MeV), the energies of initial and final ground states, respectively. Latter on we 
that the kinetic energy of the final nucleus is negligible. 

The weak interaction Hamiltonian in the neutrino mass eigenstate basis has the standard 
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form 



(u) Gf 



H^--\X) = [lL{x)^a{l+l,>^L{x)\ J^{x) + h.C. (/ = 6, /i, r), (9) 

where ja{x) is the charged hadron current. The neutrino mixing matrix Ul^^ is defined in 
Eq. (§. 

In second order of the weak-interaction we get for the (/i~,e^) conversion the following 
matrix element 



•m' 



1 



where 



with 



M 



9a 



(10) 



47r7e /-dg 

„ \ q-E^-+En-E, + iSn 

q + Ee+ + En - Ei + iSn 



(12) 



(27r)3y 2q 



n\ q-E^-+En-Ei + ien ^ 

< 0+1 E™r+a^e^^^-'-<l>(r„)|n > ■ < n| r+a^e-^^^e+^le-*^^-^^' |0+ > 



q + Ee+ + En - Ei + iSr 



(13) 



Here, TZ = tqA^/'^ is the mean nuclear radius, with tq = 1.1 fm and me is the mass of 
electron, fi is a coordinate of the ith nucleon. i?^- (^a*-) and i?e+ (/i;e+) denote energies 
(four-momenta) of the bound muon (£^^- = — ej,) and the emitted positron, respectively. 
En and are respectively energy and width of the intermediate nuclear state. $(r) is the 
radial part of bound muon in its orbit (see Appendix A). 

In the derivation of nuclear matrix element Adf.^ ^ we neglected the contribution from 
higher order terms of nucleon current (weak-magnetism, induced pseudoscalar coupling), 
which are expected to play a less important role. Following the analysis in Ref. their 



consideration can reduce the value of ■M%m„>f,e amount of about 20% by analogy to 

the Oi//5/?-decay. 
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We have normalized the nuclear matrix element ■M.%rn„>^e same way as usual 

for the corresponding Oi//5/3-decay matrix element. We note that the denominators in the 
expressions for the e"*") conversion and the Oz//3/3-decay exhibit a different behavior. This 
is because the energy of the bound muon E^- is large. The two denominators in Eq. ( jlB]) 
can be associated with the direct and the cross Feynman diagrams in Fig. [l|. One notes that 
the value of {—E^- + En — Ei) is negative. This fact implies that the widths of the nuclear 
states play an important role and that the imaginary part of the nuclear matrix element 
can be large. This point was not discussed in previous publications [0, ||, |l^, and is 
one of the motivations of our e^) conversion calculation. We want to investigate if this 
singular behavior of the amplitude can lead to an enhancement of the (/i~,e"'') conversion 
branching ratio or not. In order to simplify the numerical calculations we complete the sum 
over virtual intermediate nuclear states by closure after replacing En, Sn by some average 
values < En >, e, respectively: 



E 



\n >< n\ 



q- E^^ + En- Ei + iSn 



E 



\n >< n\ 



q + Ee+ + En - Ei + iSn 



q-E.~+ < En > -Ei + ie' 



q + Ee++ < En> -Ei + ie 



(14) 



Next we assume that the muon wave function varies very little inside the nuclear system, 
i.e., the following approximation is used 



\M 



--< 



\M 



<mv>f_ie I 



(15) 



The explicit form of < $^ >^ is given in Appendix B. 
For the width of (/i~, e+) conversion we obtain 



F 

where c 



<r^^^/>^ 



1 



TT 



E,+ K+ F{Z - 2, E,+ ) c^e < $M >' 1-^ 



nip 



(16) 



The nuclear matrix element 



= 2G'p[(mem^)/(47^m^7^)]^5f^ and ke+ = 
M 

<m^>u,e "^^^ be decomposed into the contributions coming from direct and crossed Feyn- 
mann diagrams in Fig. |l| as 



M 



(17) 



where 



dir. 



<0+| E^fcV47rE(-l)'v/2ATTjA(fce+i?fcO {rA(^^r)®n(^^i?)}oX 
kl X 

jo{qrki)jx{ke+rki/2) ^.2/2n 



n Jo q- E^- + En -Ei + ie 



[(Jk ■ o'lfA^q ) - 



9\ 



-)qdq |0j > 



kl A 



^ Mqrki)jxike+rki/2) 
n Jo q + Ee+ + En — Ei + ie 



>1 ■ o'lfliq'^ 



9l 



-)qdq |0j: >, 
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Table 1: Nuclear matrix elements of the light Majorana neutrino exchange mode of the 
(/i~,e"'") conversion in ^^Ti [see Eqs. (|17D and (0)]- The calculations have been performed 
within pn-RQRPA without and with consideration of two-nucleon short-range correlations 
(s.r.c). 







without s.r.c 






with s.r.c 
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j^cro. 

[10-2] 


[10-2] 


[10-2] 


[10-2] 


j^cro. 
[10-2] 


[10-2] 


[10-2] 


[10-2] 


0.80 


9.65 


0.23 


8.83 


13.2 


4.88 


-7.99 


4.98 


5.87 


1.00 


7.71 


3.36 


5.88 


12.5 


3.40 


-4.03 


2.37 


2.45 


1.20 


5.05 


9.09 


1.78 


14.2 


1.30 


2.71 


-1.32 


4.22 



with 

Tij Ti Tj^ T^j I T^j I , ^ j "i -^ij I -^ij I * ( ^ ) 

For the normalized nucleon form factors we use the conventional dipole form fviff) = 
1/(1 + q^klf [kl = 0.71 (GeVn fAiq') = 1/(1 + qVk\)' [Aa = 1-09 GeV]. 



4 Results and discussions 

The nuclear matrix elements of the (/x~, e"*") conversion process have been calculated within 



the proton-neutron renormalized Quasiparticle Random Phase Approximation [|18|, ^ |36 



3^. The considered single-particle model space both for protons and neutrons have been 
as follows: The full — 3hu shells plus 2si/2, Ogj/2 and Ogg/2 levels. The single particle 
energies were obtained by using a Coulomb-corrected Woods-Saxon potential. Two-body 
G-matrix elements were calculated from the Bonn one-boson exchange potential within the 
Brueckner theory. The pairing interactions have been adjusted to fit the empirical pairing 
gaps [Q. The particle-particle and particle-hole channels of the G-matrix interaction of 
the nuclear Hamiltonian H are renormalized by introducing the parameters Qpp and Qph, 
respectively. The calculations have been performed for = 1.0 and gph = 0.8, 1.0, 1.2. 



The two-nucleon correlation effect has been considered in the same way as in Ref. [^, ^ 

In calculation of the {fi~,e^) conversion nuclear matrix elements we have used the fact 
that the widths of the low lying nuclear states are negligible in comparison with their energies. 
Therefore we have carried out the calculation in the limit e ^ using the formula 

^ V--i7c6{a), (20) 



a + z£ a 

which allows one to separate the real and imaginary parts of the {fi~ , e^) conversion ampli- 
tude. 

In Table |I] nuclear matrix elements of the light Majorana neutrino exchange mechanism 
of the {fj,~,e~^) conversion in '^^Ti are presented. The adopted value of < > —Ei was 10 
MeV. We have found that our results depend weakly on this average value of the nuclear 
states within the interval 2 MeV < {< En > —Ei) < 15 MeV . However, they depend 
significantly on the details of nuclear model, in particular, on the renormalization of the 
particle-particle channel of the nuclear Hamiltonian, and on the two-nucleon short-range 
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correlation effect (s.r.c). A new feature of this (/i~,e+) conversion calculation is that the 
imaginary part of A^<m^> ^ is significant, i.e., can not be neglected. This fact was not 
noticed in the previous e+) calculations |]^, ||, |10|, pT |. 



In the further analysis we shall consider the nuclear matrix element |A1<m^>M<=| obtained 
for Qpp = 1.0 by considering the two-nucleon short range correlations. It is interesting to 
compare its value with the value of Oz//3/3-decay matrix elements for A=48 nuclear system. 
We have 



2.45 X 10 



-2 



0.82. 



(21) 



We see that the matrix element for the e^) conversion is strongly suppressed in compar- 
ison with Oi//3/?-decay matrix element by about factor of 400. It is mostly due to the large 
momentum of the outgoing positron in the (/i~,e^) conversion process. 

One can compare also the width of e"*") conversion in ^^Ti with the width of 0^(3(3- 
decay of ^^Ca. We get 



<»ni/>ee 



1.97 X 10^ 



2, Ee+) Cfj,e 



\M 



\M 



<mv>e 



\M 



\M 



< m,, > 



176. 



< > 



fie 



f22) 



The width of the {fi~, e^) conversion is enhanced mostly due to the larger available energy for 
this process. A comparison with the width of Oz//5/5-decay show that it is disfavored by smaller 
coulombic factor F(Z,E) (~ 0.623/1.8) and by significantly smaller value of associated nuclear 
matrix element (~ (0.0245/0.82)^). If we assume the effective neutrino masses < mj, >^e 
and < nil, >ee to be comparable [see Eq. (0)] we find that T^rn„> is enhanced by a factor 
of about 200 in comparison with r<m,^>^^. We have used Gqi = 8.031 x 10"^'' year~^ ||39 . 

From the experimental point of view it is interesting to compare the (/i~, e"*") conversion 
width with the width of ordinary muon capture rate. We have 



\M 



Gl [G^y + 3Gl + Gl- 2GAGp]Zf{Z, A) 



m. 



2.24 X 10"^^|A^< 



1.34 X 10"^^ 



< nil, >^e 



rrip 



(23) 



If we use the prediction for < mj, >^e coming from neutrino oscillation phenomenology, i.e., 
< >iie < 0.97 eV and more conservative bound < mj, >^e < 9 eV [see Eq. (^], we 
end up with 



4.8 X 10"^^ 4.5 X 10"^^ 



(24) 



This value is about ten orders of magnitude smaller as the estimated one for the (/i , e"*") 
conversion in ^^5* by Doi et al. There could be various reasons for this difference. First, 
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the < rriiy nuclear matrix element calculated by Doi et al. contains contributions from 
all final nuclear states and not only from the 0^^ — > 0^^ transition as in our case. It can be 
that the experimentally interesting g.s. g.s. transition exhausts only a small part from 
all allowed nuclear transitions. Second, in the simplified calculation of Doi et al. nuclear 
matrix elements were evaluated by summing the final nuclear states with closure. It usually 
leads to overestimation of the results as we know from calculation of ordinary muon capture. 
Third, the nuclear matrix element of Ref. have been evaluated by using the long-wave 
approximation. Our comparison of the (/i~,e"'") conversion and the Oi//3/3-decay (long-wave 
approximation is used) matrix elements shows that it can lead to overestimation of A^<m^>^e 
by factor up to 10^. Fourth, the problem of the ground and short-range correlations have 
been not addressed in Ref. 0]. We have found that A1<m^> ^+ matrix element for A=48 
nuclear system is strongly suppressed by both of them. It is not clear whether this effect 
is due to the chosen target. However, we note that A4<m.^> + for the conversion in ^^Ti 
consists of transition to the doubly closed shell nucleus '^^Ca, which, e.g., in the case of 
Oi//3/3-decay are less favored. To clarify this issue, the calculations of the e"'') conversion 
for other nuclear systems are necessary. 

It is worthwhile to notice that our result is in relative good agreement with the calcu- 
lations performed by Leontaris and Vergados for the (/i~,e^) conversion in ^^Ni. By 
using the same value for effective Majorana neutrino mass < rrii, >^e as in this article the 
result of Ref. corresponds to a branching ratio equal to 3.2 x 10^^^ relative to the total 
absorption of the muon for the ground state to ground state transition. The difference of 
about one order with our result for A=48 can be attributed to the nuclear physics aspect of 
the {fj;~,e~^) conversion, i.e., to a given nuclear system and the chosen nuclear model. We 



also remark that in Ref. [0 the imaginary part of the (/i , e^) conversion amplitude was 
not considered. 



5 Summary and outlook 

In summary, the lepton number violating process of e"*") conversion in nuclei have been 
studied. The light Majorana neutrino-exchange mechanism of this process have been con- 
sidered. A detailed analysis of this mode of the (/i~,e+) conversion have been performed. 
The first realistic calculation of the 0^^ 0^^ channel of this process, which is most favored 
for experimentally studies, are presented. The relevant matrix elements for A=48 nuclear 
system have been calculated within the pn-RQRPA. The effects of the ground state and 
two-nucleon short-range correlations have been analyzed. It was found that by inclusion of 
them the value of (//~,e"'") conversion matrix elements is strongly suppressed. We are the 
first, to our knowledge, showing that the imaginary part of the nuclear matrix element is 
large and should be taken into account. Further, a comparison of different relevant aspects 
with the Oz//3/3-decay process are presented. It is shown that the width of e+) conversion 
is about by factor of 200 larger as that of Oi//5/?-decay by assuming predictions for effective 
neutrino masses coming from neutrino oscillation phenomenology. Nevertheless, the studied 
neutrino exchange mode of lepton number violating {fi~,e~^) conversion is not suitable for 
experimental study being extremely small compared to the ordinary muon capture. This 
fact, however, does not disfavor further experimental study of the {fi~,e~^) conversion in 
nuclei as some other lepton number violating mechanisms, e.g., those coming from GUT's 
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and SUSY models, can dominate this process. Therefore, they should be carefully examine 
too. 
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6 Appendix A 

The bound muon wave function (IS-state) is 



where 




with = 47r/(m^e^) {a^/ae ~ rrie/m^ 5 x 10 ^) is the reduced mass of muon nucleus 
system. 



7 Appendix B 

The width for the ordinary muon capture reaction /i~ -|- (Z, A) + {Z — 1, A) can be 



written in the Primakoff form [^| 



Tm = ^K^Gf cos e,)^ < >2 Z[Gl + 3Gi + Gl- 2GAGp]f{Z, A), (27) 
where muon average probability density over the nucleus is 

< ^ ' 7 ' — . 28 

j p[x)d-^x 

p{x) is the nuclear density. To a good approximation it has been found 

< >^= (29) 

i.e., the deviation from the behavior of the wave function at the origin has been taken into 
account by the effective proton number Z^ff. The values of this effective charge has been 
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calculated for the nuclear systems of interest in Ref. [1^. In particular, one finds Z^ff = 17.6 
for Z = 22. The quadratic combination of the weak coupling constants is 

[G^ + 3G\ + Gl- 2GaGp] ^ 5.9. (30) 

The function f{Z,A) takes into account the two-nucleon correlations given by ||lT| 



liA, Z) = 1 - O.O3A + o.25( A _ 1) + 3.24(A - i - I ^ - (31) 
This Pauli blocking factor for ^^Ti takes the value /(22,48) = 0.11. 
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